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Abstract 

We study dynamical behaviors of the dark energy models interacting with dark matter and un- 
particle in the standard flat FRW cosmology. We considered four different interacting models and 
examined the stability of the critical points. We find that there exist late-time scaling attractors 
corresponding to an accelerating Universe and the alleviation of the coincidence problem depends 
on the choice of parameters in the models. 
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I. INTRODUCTION 



Many observations confirm that we live in an accelerating Universe. Within the framework of Einstein 
gravity, this acceleration can be explained by the mysterious energy component named dark energy (DE) with 
negative pressure, which occupies almost 70% of the content of the universe at present 3,3,0, 4 1. The simplest 
explanation for dark energy is the cosmological constant, which is a term that can be added to Einstein's 
equations. This term acts like a perfect fluid with an equation of state (EoS) lu x = —1, and the energy density 
is associated with quantum vacuum. Although this interpretation is consistent with observational data, at 
the fundamental level it fails to be convincing. The vacuum energy density is far below the value predicted 
by any sensible quantum field theory, and it suffers the coincidence problem, namely, why the DE and the 
Dark Matter (DM) are comparable in size exactly right now. To overcome the coincidence problem, some 
sophisticated dynamical scalar field DE models, such as quintessence, k-essence and phantom field, have been 
put forth to replace the cosmological constant [5|. However, these scalar field DE models can not still resolve 
the coincidence problem. 

The interacting DE model has been regarded as a possible way to alleviate this problem [f| • In these models 
it is assumed that there exists a nonzero interaction between DE and DM in the Universe that gauges DE 
transfers to DM which allows us to create an equilibrium balance in the evolution of the Universe, so that 
the density of DE keeps the same order as that of DM at late times. Although lacking of a fundamental 
explanation for the source of the interaction, the interacting models arc manifested as the useful and robust 
models at the order of one standard deviation by some data sets from observational cosmology including 
( MB shift parameter. BAO, age parameter and supernova observations and so on [6|, \[, 8L Therefore, the 
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interac ting models have attracted a great deal of interest 
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24, 



25j . Recently, some attempts have been 



from a fundamental field theory point of view 



proposed to describe the interaction between DE and DM 
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The investigations above have been focused only on the interaction between DE and DM. However, it is 
physically reasonable and even expected from a theoretical point of view, that DM as well as DE can interact 
with other dark components of the universe 
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29(. Recent investigations shows that unparticle can be 



treated theoretically as another important dark components in the Universe. The concept of unparticle is 
introduced first by Georgi 3CJ , which is based on the hypothesis that there could be an exact scale invariant 



hidden sector resisted at a high energy scale (for a recent review of unparticles, see 



33 |34j). Although the 



fundamental energy scale of such a sector is far beyond the reach of todays or near future accelerators, it is 
possible that this new sector affects low energy phenomenology. These effects are described as unparticle in 
the effective low energy field theory because that the behaviors of these new degrees of freedom are quite a 
different those of the ordinary particles. For example, their scaling dimension does not have to be an integer 
or half an integer. Recently, a lot of work 
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32j have been focused on the new collider signals for unparticle 
physics. One of interesting feature of unparticle is that it does not have a definite mass and instead has a 
continuous spectral density as a consequence of scale invariancc 



p{P 2 ) = A d J(P°)8(P 2 )(P 2 ) 
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where P is the 4-momentum, Ad u is the normalization factor and d u is the scaling dimension. The theoretical 
bounds of the scaling dimension d u are 1 < d u < 2 (for boson unparticle) or 3/2 < d u < 5/2 (for fermion 



unparticle) 



34j . The pressure and energy density of the thermal boson unparticle are given by 
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where C(d u ) = B(3/2,d u )T(2du + 2)((2du + 2), while B, T, £ are the Beta, Gamma and Zeta functions, 
respectively. Thus, the EoS of boson unparticle reads [35 1 

For the fermion unparticle, we find the EoS has the same form as that of boson one. Obviously, the EoS of 
unparticle ui u is positive which is different from that of DE and DM. Since the unparticle interacts weakly with 
standard model particles, it can be regarded as a new form of dark component. Recent inve stig ations show 
that the unparticles play an important role in the early universe 



36] and black hole physics 



3J. Therefore, 



it is natural to ask whether there exists some new properties in the late-time evolution of the Universe if the 
unparticle takes part in the interaction with DE and DM? In this paper we choose four different coupling 
terms and study the dynamical behaviors of the interacting DE models with DM and unparticle by the phase 
space analysis method to discuss further stability of the critical points and their cosmological implications. 

The paper is organized as follows: in sections II, we construct a cosmological scenario in which DE interacts 
with DM and unparticle, and then we present the formalism for its transformation into an autonomous 
dynamical system which is suitable for a phase space stability analysis. In Sec. Ill, we consider some special 



coupling forms and perform the phase space analysis of the corresponding interacting DE model with DM and 
unparticle, and then discuss their cosmological implications. Our conclusions and discussions will be presented 
in the last section. 

II. INTERACTING DARK ENERGY MODEL WITH DARK MATTER AND UNPARTICLE 

In the Einstein theory, the a flat FRW universe is described by the standard Friedmann equation and 
Raychaudhuri field equation 

H 2 = £p, (4) 
H = ~{p + P). (5) 

H is the Hubble parameter and k is the constant 8irG. The total energy density p = p m + p x + p u , where p m , 
p x and p u correspond to the energy densities of DM, DE and unparticle, respectively. For simplicity here we 
have neglected the radiation and baryons since we are concentrating on the late time accelerating Universe. 
The interaction among DE, DM and unparticle can be described in the background by the balance equations 

p x + 3H(l + u x )p x = Ti, 

Pm + 3Hp m = r 2 , 

Pu + 3H(1 + uj u )p u = T 3 (6) 

Here the terms T\, T 2 and T 3 describe the coupling among DE, DM and unparticles. The total conservation 
equation demands that 

ri + r 2 + r 3 = o. (7) 

To analyze the evolution of the dynamical system, we introduce the dimensionless variables 

= Kp x _ Kp m _ KPu d ,m 

X ~3H 2 ' y ~3H 2 ' Z ~3H 2 ' dN~Hdt' [) 

where N = In a is the number of e-folding to represent the cosmological time. Using the above definitions, the 

Hubble equations can be rewritten as 

x + y + z = l, (9) 

and 

UJ X X + UJ U Z 
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1 + 



x + y + z 
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= --(1+u) x x + uj u z). (10) 



The effective total EOS ujtot is given by 

LU x p x +LU u p u LU X X + LU U Z 

u to t= ; ; = ; ; = lu x x + lu u z. (11) 

Px + Pm + Pu X + y + Z 

Once the concrete forms of the coupling terms Y\, T% and 1^ are given, the equations of motion ((4]), ([5]) and 
([5]) can be transformed to an autonomous system containing the variables x and y and their derivatives with 
respect to -/V = In a. This autonomous system has generally the form: X' — f(X), where X is the column 
vector constituted by the auxiliary variables, f(X) is the corresponding column vector of the autonomous 
equations, and prime denotes derivative with respect to N — In a. Solving the equations X' — 0, we can 
obtain the critical points X c . Then, in order to describe the stability properties of these critical points, we 
can expand the equations X' — f(X) around X c by setting X = X c + SU with SU the perturbations of 
the variables considered as a column vector. Thus, for each critical point we can expand the equations for 
the perturbations up to the first order as: 5U' = 3 • SU , where the matrix S contains the coefficients of the 
perturbation equations. Through the analysis of the eigenvalues of S, we can describe the stability of each 
critical point. In general, for an arbitrary coupling terms T; it is difficult to obtain the analytical forms of the 
critical points. Here we only consider some specific forms of I\ and examine the stability of interacting DE 
models with DM and unparticlc in the next section. 

III. ANALYSIS OF STABILITY IN THE PHASE SPACE 

In this section, we will consider some concrete coupling forms Tj and analysis the stability of the corre- 
sponding interacting DE models with DM and unparticle, and then discuss their cosmological implications. 

A. Interacting model I 

In general, there is as yet no basis in fundamental theory for a specific coupling in the dark sectors because 
that the nature of dark sectors remain unknown. Thus, all coupling models discussed at the present moment 
are necessarily phenomenological [6!]. There are two criterions to determine whether some models can be 
more physical justification than the others. One is to confront observations. The other is to examine whether 
the coupling can lead to accelerated scaling attractor solutions [15j, which is a decisive way to achieve similar 
energy densities in dark sectors and alleviate the coincidence problem. Motivated by analogy with dissipation 
of cosmological fluids, the coupling terms T, which are proportioned to the densities p x and p m have been 



studied in the context of quintessence [6j and phantom 
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221 models. In this section, we first 



consider the case in which all of the coupling terms I\ are proportioned to the density of DE p x . Here we 
choose T 2 = T 3 so that we can obtain the analytical expression for the critical point, which is very convenient 
for us to study the dynamics of the model in the following calculations. The concrete expressions of I\ are 

ri = -6bH Px , r 2 = r 3 - 3bH Px , (12) 

where 6 is a positive coupling constant. It means that in this model only DE can convert to DM and unparticle. 
Using dimensionless variables, the dynamical equations of the system can be expressed as 

x' = —6bx — 3ui x x + 3x(u x x + u u z), 
y' = 3bx + 3y(w x x + u> u z), 

z' = 3bx — 3oj u z + 3z(w x x + ui u z). (13) 

Solving the equations x' = 0, y' = and z' = 0, we obtain the critical points (x c , y c , z c ): 

■ Point Ai : (0, 1, 0), 
• Point B x : (0, 0, 1), 

_ p oint c . f {2b + u x ){2b + lo x - u u ) -b{2b + u; x -Lu u ) -b{2b + u x ) \ 

\ (u x - u u ){b + uj x ) + buj x ' (w x - ui u )(b + Wj) + bw x ' (w x - iv u )(b + co x ) + bu x J ' 

The points Ai and B\ imply that our Universe is dominated by DM and unparticle, respectively. The condition 
for the existence of the point C\ is < b < — The critical point C\ denotes that the DE, DM and unparticle 
can be coexisted in the late-times of the Universe. 

After some operations, we can obtain the 3x3 matrix E of the linearized perturbation equations 



u> x (6x c — 3) + 3u u z c — 66 3lo u x c 

3(u> x y c + b) 3(lu x x c + lu u z c ) 3uo u y c 

3{u) x z c + b) 3uj x x c + u> u (6z c - 3) 



(15) 



The eigenvalues of the coefficient matrix S encode the behavior of the dynamical system near the critical 
points. If the real parts of all of the eigenvalues of the matrix El are negative, then the critical point is a stable 
point, otherwise it is unstable. Through some careful calculations, we obtain the eigenvalues of the coefficient 
matrix for these critical points 

• Point Ax : X 1 = 0, A 2 = -3w u , A 3 = -3(26 + lu x ), 

■ Point B x : \ x = 3lo u , A 2 = 3uj u , A 3 = -3(26 + u) x - us u ), 

■ Point Ci : Ai =3{2b + u x ), A 2 = 3(26 + lu x ), A 3 = 3(26 + uj x - w„). (16) 
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For the point Ai, the eigenvalue Ai is non-negative, which indicates that Ai is not a stable point. For the point 
Bi, we also find that both of the eigenvalues Ai and A2 are positive because that the EoS of the unparticle 
u>u > 0. Therefore the point B\ is not the stable point. For the point Ci, when < 6 < — all of the 
eigenvalues Ai, A2 and A3 are negative, which means that C\ is a stable point as shown in figure (1). The 
stable region of C\ is not affected by the EoS of the unparticle, but the position of C\ in the phase space is 
decided together by u> x , ui u and the coupling constant 6. From Eq. lfTTj) . we also learn that the effective total 
EoS at point C\ is Lu to t = 26 + w K , which is independent of lu u . When b —> —u> x /2, we find that uj to t tends to 
zero. It is consistent with the result shown in figure ([2]) in which as the coupling is strong enough the Universe 
will be dominated by DM. Moreover, for the critical point C%, we find that in its stable region none of the 
coordinates x c , y c and z c in phase space vanishes, which means that the coincidence problem can be alleviated 
in the Universe described by the point C\ . Since the effective total EoS at point C\ is uJtot = 2b + uj x we else 
obtain a oc — (3w x + 66+1) 1 h^+^+i) _2 anc [ p ^ a -3(w x +2b+i) ^ means that point C\ is an accelerated 
scaling solution as 6 < — ^ — | and there is singularity in the finite future as 6 < ^ — -|. 




FIG. 1: The phase diagram of interacting dark energy with DM and unparticle through the coupling terms (|12[) . The 
point Ci is the critical point. Here we choose the values uj x = —1.2, uj u = 0.28 and b — 0.5 in the stable region 
b < -w x /2. 

B. Interacting model II 

In the previous discussions, we studied only the case DE converts to DM and unparticle and do not consider 
the exchange of energy between DM and unparticle. In this section, we will consider another special case in 
which besides DE transfers to DM, DM can also be converted to unparticle. Similarly, in order to obtain the 
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FIG. 2: Variety of x c , y c and z c with b at the critical point C\ for fixed ui x = — 2 and uj u = 0.3. The coupling constant 
b is located in the stable region < b < —u x /2. 



analytical forms for the critical points, we assume that the coupling terms I\ have the forms 

Ti = -3bHp x , T 2 = 3bH(p x - p m ), T 3 = 3bHp m , 

respectively. 

The dynamics of the system can be described by 

x' = — 3bx — 3uj x x + 3x(lj x x + u) u z), 
y' = 3b(x-y) + 3y(ui x x + uj u z), 
z = 3by — 3lo u z + 3z(u x x + u u z), 
and the critical points (x c , y c , z c ) are 

• Point A 2 : (0, 0, 1), 

• Point B 2 : (0, I-—, — ), 

0J U U! u 



Point C 2 : 



u x (b + ui x - iv u ) -b{b + u x - uj u ) 



(17) 



(18) 



(19) 



y uul+u u {b- uu x y uj% + ui u (b - u x ) ' ul+ uo u {b - uu x ) )' 
The condition that the point Ci exists is < b < uo u — ui x . Similarly, the 3x3 matrix S of the linearized 



perturbation equations reads 

-3(lu x + b - 2uj x x c - lo u z c ) 3lu u x c 

3u> x (y c + b) 3(uj x x c + lu u z c ) - 3b 3uj u y c 

3u> x z c 3b 3lu x x c + lu u (6z c - 3) - 

The eigenvalues of the coefficient matrix for these critical points are 

• Point A 2 : Xi = 3uj Ui X 2 = 3(u> u -b), A 3 = -3(6 + lo x - u) u ), 

■ Point B 2 : \\ = 36, \ 2 = —3uj x , A 3 = 3(6 — u u ), 

■ Point C 2 : Ai = 3u> x , \ 2 = 3(u> x + 6), A 3 = 3(6 + lo x - u u ). 



(20) 



(21) 
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For point A2, since u) u > 0, the eigenvalue Ai always is positive, which indicates that A2 is an unstable point. 

Moreover, the EoS of DE cu x < and the coupling constant b > means that both of the eigenvalues Ai and A2 

of the point B2 are positive. Thus the point B2 is an unstable point. For the point C2,when — uj x < b < lu u —u x , 

the sign of A2 is always opposite to the signs of Ai and A2, which leads C2 to a saddle point. However, when 

< b < — tu x , we find that all of eigenvalues (Ai, A2 and A3) are negative, which indicates that C2 is a 

stable point as shown in figure (3). Meanwhile, we also find none of the coordinate components of C2 in the 

phase space disappears, which means that in this case there exist three components (DE, DM and unparticlc) 

in the late-time universe and the coincidence problem can be alleviated. From Eq. (|ll[) , we learn that the 

2 9 

effective total EoS at point C2 is Lo to t = b + u x . Therefore we obtain d oc —(1 + 3b + 3lo x ) t a v-+ b + u *'> and 
p oc a - 3 ( 1 + fc + w *) , it means that point C2 is an accelerated scaling solution as b < —lu x — 1/3 and there exists 
singularity in the finite future as b < —oj x — 1. For fixed uj x , the effective total EoS increases with b. Figure (j4|) 
shows that in the Universe described by C2 the density of the unparticle increases and DE decreases with the 
coupling constant b, while the density of DM first increases and then decreases. This means that the Universe 
will be dominated by unparticle if the coupling constant b is large enough. 



FIG. 3: The phase diagram of interacting dark energy with DM and unparticle through the coupling terms (I17p . The 
point Ci is the critical point. Here we choose the values lo x — —1.2, u u = 0.28 and b — 0.5 in the stable region b < — u x . 
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C. Interacting model III 



In Ref. 2l| , we have investigated the dynamics of the interacting DE models in which the coupling terms Tj 
contains the product of the densities of DE and DM and find that this new type of dark sector coupling leads 
to the more interesting accelerated scaling solutions and presents us more complicated features in dynamical 
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phase space. Therefore, in the following sections, we will consider the cases in which the coupling terms 
Ti contain the product of the density of DE, DM and unparticle and to see whether it presents some new 
properties or not in the evolution of the Universe with the unparticle component. For mathematical simplicity, 
we assume firstly the interaction among DE, DM and unparticle have the forms 

Ti = -6bnH~ 1 p x p u , r 2 = T 3 = 3bnH~ x p x p u . (22) 

As in the model (fT2"]) . these coupling terms also denote that DE can be transfer to DM and unparticle. The 
dynamical equations of the system can be written as 

x' = — Qbxz — 3uj x x + 3x(u> x x + u> u z), 

y' = 3bxz + 3y(uj x x + lj u z), 

z' = 3bxz — 3lu u z + 3z(lu x x + u> u z), (23) 
and the critical points (x c , y c , z c ) arc 

• Point A 3 : (1, 0, 0), 

• Point B 3 : (0, 0, 1), 

Point C ■ ( UJu ^ b + UJx -^u) {lo x - u u ) 2 + 3b{uj x - u u ) + 2b 2 b(b + uj x - uj u ) \ 
° m 3 ' \b (2b + 2uj x -Uu)' b (2b + 2oj x - uj u ) ' b (2b + 2uj x - lo u ) ) ' ^ ' 

The condition that the point C3 exists is b > uj u — lo x . Repeating the previous operations, we can obtain the 
3x3 matrix 5 of the linearized perturbation equations 



uj x (§x c — 3) + 3uj u z c — 66 3lu u x c — 6b 

3(u> x y c + bz c ) 3(uj x x c + lo u z c ) 3(uj u y c + bx c ) 

3(lu x + b)z c 3to x x c + to u (Qz c — 3) + 3b 



(25) 
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and the eigenvalues of the coefficient matrix for these critical points 



Point Ai 



Ai = 3u x , A 2 — 3uj x , A 3 = 3(6 + uo x - lo u ), 



• Point B 3 : Ai = 3uj u , A 2 = 3u u , A 3 = -3(26 + u x - uj u ) 7 

3uj u u) x 



Point C 3 : Ai 



26 + 2u) x — uj u 



A 2 ,3 — 



3(buj u u! x ± y/ (buj u uj x [16b 3 + 8b 2 (5m x - 4u u ) + A(uj x - uj u ) 2 (2uj x - u u ) + 6(32ug - hlu x uj u + 20cjg)]) 



26(26 + 2uj x - uj u ) 



(26) 



Through the similar analysis, we find that A3 is a stable point when 6 < oj u — lo x and is a saddle point 




FIG. 5: The phase diagram of interacting dark energy with DM and unparticle through the coupling terms (|22[) . In 
the left figure, we fix the values (u x — —1.2, uj u = 0.28 and b = 0.5) which meets b < u) u — u x and the critical point is 
A3. In the right, we choose the values (u x = —1.2, u u — 0.28 and 6 = 1.8) which satisfies b > u> u — ui x and the critical 
point is C3. 



when 6 > lu u — u x . The point B3 is an unstable point because that lo u > 0, which leads to the signs of the 
eigenvalues Ai ;2 always are positive. For the point C3, it is a stable point when 6 > lu x — lu u . We can obtain 
the total effective EoS u>tot — 26+2^"-^ • ^ depends not only on the EoS of DE and unparticle, but also on 
the coupling constant 6. From the total effective EoS uJtot, it is easy to obtain that point C3 is an accelerated 
scaling solution as oj u — to x < 6 < (u> u — 2lu x — 3uj x lo u ) /2. While for 6 > (cu u — 2uj x — 3u) x u) u ) /2, we find 
that Wt t > —1/3, which means that it is a decelerated scaling solution. This can be explained by that the 
larger 6 denote the more DE transfers to DM and unparticle in the Universe. While, as 6 < lo x — ui u , the 
point C3 is an unstable point. Moreover, in this model, when 6 is small the coupling the Universe will enter 
the era dominated by DE because A3 is a stable point (as shown in the left figure in Fig. (5)). When 6 is 
larger the Universe will enter a stable stage described by point C3 (as shown in the right figure in Fig. (5)), 
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which contains the DE, DM and unparticle. Thus, when b > lu u — lo x this model can resolve the coincidence 
problem. Moreover, from the Fig.©, we also find with the increase of b the density of DM increases in the 
Universe, but both of DE and unparticle decrease. This implies that the Universe will be dominated by DM 
when the interaction is very strong. Figure also tells us that with the increase of the coupling constant b 
the EoS Lo t ot increases so that the Universe can make a transition from an accelerating expansion phase to a 
decelerating one. As b tends to infinity the EoS Lo to t approaches to zero, which is consistent with that of DM. 
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0.6 
0.4 
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2 3 4 5 6 7 8 

b 

FIG. 6: Variety of x c , y c and z c with b at the critical point Cz for fixed u> x — —1.7 and u) u = 0.3. Here the coupling 
constant b is located in the region b > lu u — ui x . 
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FIG. 7: Variety of the effective EoS uitot with b at the critical point C3 for different ui x and ui u . Here the coupling 
constant b is limited in the region b > ui u — ui x . 

D. Interacting model IV 

In the model (III), we consider only the case the coupling terms do not contain the density of DM. In the 
following model, we will consider the interaction includes p m and suppose that Ti have the forms 

Ti = -3bKH~ 1 p x p u , T 2 = ibnH^ 1 (p x p u - p m p u ), T 3 = 3bnH^ 1 p m p u . (27) 
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Although the coupling forms are more complicated than the model it also presents us the analytical 

expressions for the critical points which brings us the simplicity in the later calculation. This type of Tj also 
denote DE can be transfer to DM and at the same time DM can also be convert to unparticle. The dynamical 
equations of the system are given by 

x = — 3bxz — 3u> x x + 3x(lu x x + lu u z), 
y = 3b(xz - yz) + 3y{uj x x + uj u z), 

z = 3byz — 3u> u z + 3z(u> x x + u> u z), (28) 
and the critical points (x c , y c , z c ) can be expressed as 

• Point A 4 : (1, 0, 0), 

■ Point B 4 : (0, 0, 1), 

„ . , „ fb + u x -u> u (u> u -uj x )(b + u x ~ u u ) u x (uj x ~u) u )\ 

■ Pomt G 4 : , — — r , — — — . (2\)) 

\ b b (b - LU U ) b (b- LU U ) J 

The critical point C4 exists only when b > lo u — lu x . Similarly, we can obtain the 3x3 matrix S 



u x (6x c — 3) + 3lo u z c — 3b 3{uj u — b)x c 

3(u x y c + bz c ) 3(lu x x c + lu u z c - bz c ) 3(u u y c + bx c ~ by c ) 

3uj x z c 3bz c 3uj x x c + w„(6z c - 3) + 3by c 

The eigenvalues of the coefficient matrix for these critical points are 

• Point A 4 : Ai = 3lo x , A2 = 30^, A3 = 3(uj x — u> u ), 

■ Point B 4 : Ai = 3u> u , A 2 = 3(uj u - b), A 3 = -3(6 + uj x - uj u ), 



(30) 



3uj x (b + lu x - 2lu u ) 3(bu x ± ^buj x [A{uj x - uj u ) 2 + b(5uj x - 4w„)]) 

• romt O4 : Ai = : , A2.3 = tt, • (oi) 

b — uj u Lb 

Obviously, all of eigenvalues (Ai, A2 and A3) for the point A4 are negative because the EoS of DE io x < 0, 
which indicates that in this case the point A4 always is a stable point. The point B4 is an unstable point 
since Ai is positive. From the previous discussion we find that the point C4 exists only when b > lo u — lo Xi 
which leads to that A2 > and A 3 < 0. This means that the point C4 is a saddle point. Since the point A4 
describes a Universe filled with only DE, the interaction (|27|) cannot resolve the coincidence problem. 



IV. CONCLUSIONS AND DISCUSSIONS 



In this paper we have studied the dynamical behaviors when DE is coupling to DM and unparticle in the 
standard flat FRW cosmology. We considered four different interacting models and examined the stability of 
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critical points. In all the examined models we found that there exist a stable late-time scaling solution which 
corresponds to an accelerating universe. This feature was expected since DE cosmology has been constructed 
in order to always satisfy this condition. Moreover, for all the studied cases, we also find the accelerating 
stable solutions if we choose the appropriated parameters in the models. Except the fourth model, the stable 
solutions in other models admit that DE coexists with DM and unparticlc in the Universe. Our result also 
implies that in these models the coincidence problem can be alleviated only when the coupling is not strong 
enough. 
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